Crossovers in ScalePree Networks on Geographical Space 
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Complex networks are characterized by several topological properties: degree distribution, clus- 
tering coefficient, average shortest path length, etc. Using a simple model to generate scale-free 
networks embedded on geographical space, we analyze the relationship between topological prop- 
erties of the network and attributes (fitness and location) of the vertices in the network. We find 
there are two crossovers for varying the scaling exponent of the fitness distribution. 

PACS numbers: 89.75.Hc; 89.75.Da 



Many natural, social, and technological systems can be 
described in term of complex networks, in which vertices 
represent interacting units, and edges stand for interac- 
tions among them 0, S 0- The complex networks, 
which are far from absolutely regular or complete ran- 
dom, are characterized by several properties: degree dis- 
tribution, clustering coefficient C, average shortest path 
length L, etc. Many real networks exhibit a scale- free 
degree distribution P(k) ~ k~" , typically with scaling 
exponent 2<^<3[l[|l|l|l|l|i|. Most real networks 
have a large clustering coefficient, which is defined as the 
probability that a pair of vertices with a common neigh- 
bor are also connected to each other [jj. The local clus- 
tering coefficient usually decreases with the degree ||. 
In addition, small-world effect is seen in many networks 
0- For many network models, the average shortest path 
length grows logarithmically L tx InN or more slowly 

In order to understand the structure of the complex 
networks, many models have been proposed. Such mod- 
els can be grouped into two main classes. The first class 
contains growing networks with preferential attachment. 
Barabasi and Albert (BA) proposed this type model orig- 
inally 0. For BA model, the scaling exponent v is al- 
ways 3 and its clustering coefficient is relatively small 
for large size. Then, several modified models have been 
presented to reproduce the realistic aspects of networks 
0. 0. [HI Hill The second class contains static networks, 
where each vertex has a intrinsic fitness measuring the 
its importance or rank 0, 0, Q, 0, Q^. In several 
models, the location in geographical space is also taken 
into consideration [13 E3 Ea 123, [111, (13, (3, (2^1 . In 
such models, the topological properties of the network 
are essentially determined by the characteristics of the 
vertices. The purpose of this paper is to make clear the 
relationship between the topology of the network and the 
attributes of the vertices for a simple model of the second 
class. We find when the scaling exponent 7 of the fitness 
distribution varies, there are two crossovers at 7 = 2 and 
7 = 3. 

Our model is defined in the following. We consider 
N vertices. We assume that each vertex has an fitness 
a,i(i = 1,2,..., N). For simplicity, the fitness values are 



assigned deterministically as 



(i = l,2,...,N). 



(1) 



for 7 > 1 The case of 7 ~ 2 is known as Zipf 's law. When 
N is adequately large, the distribution of the fitness is 
given approximately as 



p(a) = (7 - l)a-T + 
in the finite support 



5(o-l) S(a-N—) 



2N 



2N 



1 < a < N~. 



(2) 



(3) 



Here 6(x) denotes Dirac's delta function. Thus, the dis- 
tribution of the fitness follows the power law p(a) oc a~ 7 
with slight adjustments at the both side. In addition, the 
vertices are distribute randomly in rf-dimcnsional space 
with uniform distribution. For simplicity, the distance 
is defined by L-max norm, and the boundary condition 
is periodic. We assume that the fitness and the location 
are independent mutually. The condition to link vertices 
i and j is 



didj 



< 



(4) 



where denotes the distance between these vertices 

and 9 is a threshold. Here, the threshold value 9 is cho- 
sen so that the total number of connections equals mN. 
Thus the average degree is given by (k) = 2m. The net- 
work resulting from our method has a scale-free degree 
distribution as is shown in Fig. 1(a). 

Since the vertices follow the uniform distribution in the 
unit d-dimensional cube, the probability to link a pair of 
vertices with fitness a and a' is given as 



r(a, a') = min(0aa', 1). 



(5) 



The average degree for a vertex with fitness a is calcu- 
lated as 



k(a) = N r(a,a')p(a')da'. 



(6) 



2 



Inserting J5J and (|5j into ©, we obtain the approximate We can estimate the threshold value 9 from the fact that 
form for large N the average degree is described as 

2(7- l)N -jN^t 



7-2 



9a (a < 9~ 1 N~ 



N (a>9- 1 N—) 



(7) 



2m = j k(a)p(a)da . (8) 
Inserting J2J and Q into ©, we get for large N 



2m ~ N 



2(7- l) 2 9 - (3 7 2 - 5 7 + 2)6W^ + ( 7 2 _ 7 )6l7-i _4( 7 _ 2)6>T- 1 lniV- 2( 7 2 -3 7 + 2)6lT- 1 ln( 

2(7-2) 2 



r 



(9) 



The asymptotical solution 9 for large N is described as 



2m f 7 
~N 



7-1, _ 
2m(2-7)\ t- 1 
AlniV 



( 7 >2) 



(1<7<2). 



(10) 



Thus, the asymptotical behavior changes at 7 = 2. For 
7 > 2, the asymptotic form of Eq. (JJJ is given as 



&(a) ~ 2m 



7 



7-1 



a (7 > 2). 



(11) 



Accordingly, fc(a) is proportional to a. On the other 
hand, for 1 < 7 < 2, Eq. J7J is approximately 



k(a) 



2m 
InN 



f- 1 (1< 7 < 2) 



(12) 



for a > 6~ 1 N : y= r . Thus, fc(a) follows a power law decay 
with exponent 7 — 1. 

Let us now calculate the degree distribution. The de- 
gree distribution is calculated as 



P(k) = I P(k\a)p(a)da. 



(13) 



Here, the conditional probability P{k\a) that vertex with 
fitness a has degree k is given by binominal form: 



P{k\a) 



k(a) 



1 - 



k{a) 



N-k 



(14) 



For 7 > 2, if the inside of the integral of (|13|l has the 
maximum in the range ©, the integral is approximated 
by using the gamma function. Accordingly, in the region 



2m 



7 



7-2 _l 

7 < k < 2m- N-1- 1 + 7, 



7 — 1 7 — 1 

the degree distribution (|13|l is described as 



(15) 



(16) 



(7- l)''- 2 k\ T(A-7 + 2)' 




10 100 1000 1000c 

degree: k 



■'■'••■•'••■•*• 

(b) 

1.5 2 2.5 3 3.5 4 
exponent of fitness distribution: 7 




10 100 1000 10000 

degree: k 



1.5 2 2.5 3 3.5 4 
exponent of fitness distribution: 7 



FIG. 1: (a) The degree distribution P(k) obtained numeri- 
cally for m = 3, d = 2, N = 10000, and 7 = 1.5 (circles), 
7 = 2.5 (squares) and 7 = 3.5 (triangles). This data is aver- 
aged over 100 realizations, and the bin is taken logarithmically 
to reduce noise. The solid curves stands for the theoretical 
prediction 1171 and 11911 . (b) The exponent values calculated 
using the maximum likelihood method are shown as a func- 
tion of 7 for m = 3, d = 2, N = 10000. (c) The average 
nearest neighbor degree (ANND) for the same parameters as 
in (a). The solid curves correspond to theoretical results 1211 1. 
122H . and 123H . (d) The assortativity r (degree correlation) ob- 
tained numerically is shown as a function of 7 for the same 
parameters as in (b). 



For k 3> 1, we get the scale- free degree distribution: 

P(k)^(2my- l{ ^—^k^ ( 7 >2). (17) 

In this case, the scaling exponent equals that of the 
fitness distribution. For 7 < 2, if the inside of this 
integral of Eq. i|13|) has the maximum in the range 



3 



(8~ 1 N^ rT , N~), the degree distribution is calculated 
in the same way to the case of 7 > 2. Accordingly, in the 
region 



27 , 2mN 
< k < -^^ + 2, 



2-7 In TV 

the degree distribution is described as 



P(k) 



2m 

TnN 



k~ 



(1<7<2) 



(18) 



(19) 



As a result, in this case, the degree distribution is inde- 
pendent of 7, and the exponent is always 2. These anal- 
yses are consistent with the numerical results (see Fig. 1 
(a) and (b)). Note the degree distribution is independent 
of the dimension d in the both cases. 

In addition to the degree distribution, we study the 
degree-degree correlation P(k'\k), which measures the 
probability of a vertex with degree k to be linked to a 
vertex with degree k! . In order to characterize this cor- 
relation, it is useful to work with the average nearest 
neighbor degree (ANND), which is defined as k nn (k) = 
Y2 k , k'P(k'\k) 26]. Before estimating k nn (k), we esti- 
mate the ANND of a vertex with fitness a, which is cal- 
culated as 



(20) 



Here the last term adding one is due to the fact that 
the nearest neighbor vertex has at least one connection. 
Eliminating a with using (|11J) . we obtain an approxima- 
tion for k nn {k). For 7 > 3, we obtain the asymptotical 
form for large N 



r(a, a')k{a')p{a')da! 

k„ n {a) = - — - r + 1 

r(a, a')p(a')da' 



k n n(k) 



2m(7 - 2) 2 
(7-l)(7-3) 



+ 1 (7 > 3). 



(21) 



Thus, the ANND k nn (k) is independent of k. This result 
indicates there is no correlation between degrees of linked 
pairs. However the numerical result shows there is a small 
positive correlation (see Fig. 1 (c)). This correlation may 
be due to the fluctuation of the vertex density in the d- 
dimensional space. For 2 < 7 < 3, the asymptotical form 
of ANND is 



knn (&) 



A 



A 



7+1 3-7 

-N~ - 1 



1 



2(7-1) 

(fc< 2+1^,2 < 7 < 3) 

(k> 2^A^,2< 7 <3) 



where A 



7-2 

(22) 

2m( 7 -2) 2 _ ( 7 -l) 3 ' T nnH a _ 7(3-7) 

( 7 -l)(3-7)' a ~ ( 7 -2) 4 --' ' dUU P — 2(7-2)^ ' 

This result indicates that the ANND is constant for small 




d = I 

d = 2 
d = 3 



(b) 



degree: k 



d= 1 
d = 2 
d = 3 



(d) 



degree: k 



exponent of fitness distribution: y 



FIG. 2: The clustering coefficient as a function of k for m — 3, 
N = 10000, and 7 = 1.5 (a), 7 = 2.5 (b), and 7 = 3.5 
(c), where circles, squares and triangles correspond to the 
numerical results for d = 1, d = 2, and d — 3, respectively. 
In (a), these three plots coincide exactly. The solid curves 
correspond to the predictions from the numerical integrate of 
12 11 . In (d), the average clustering coefficient as a function of 
7 for m = 3, N = 10000 and d = 1 (circles), d = 2 (squares), 
and d = 3 (triagles) . 



k and decays approximately k nn (k) cx k ( 3 7 ^ for large 
k. For 1 < 7 < 2, we obtain 



1 



k n n{k) — < 



C 2miV(2-7)(27-l) 
7 In 

(k < 1/(2- 7), 1< 7 < 2) 

4mA ^ fc + ln(2-7) + 7-l/2 + ! 
(2fc-l)lnAT 
(fc > 1/(2- 7), 1< 7 < 2) 



(23) 

This result indicates that the ANND decays k nn (k) oc ^ 
for large k. Figure 1 (c) shows these analyses agree well 
with the numerical results. In addition, we calculate nu- 
merically the assortativity r defined by Newman |27l ] for 
several values of 7. The assortativity denotes Pearson 
correlation coefficient of the degrees at either ends of an 
edge. Figure 1 (d) shows the network has positive (neg- 
ative) degree correlation if 7 > 3 (7 < 3). 

The clustering coefficient is calculated as a function of 
a as follows 



C(a) 



r 3 (a,a',a")p(a')p(a")da / da" 



r(a, a')r{a, a")p{a')p{a")da'da" 



(24) 



where r 3 (a, a', a") denotes the probability that three ver- 
tices with fitness values a, a', and a" form a triad. Since 
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FIG. 3: (a) Log-linear plot of the average shortest path length 
L vs. the number of vertices N for different values of 7 (7 = 
4.0,3.5,3.0,2.5,2.0,1.5 from top to bottom) and d = 2. (b) 
The same data in log-log plot. 



it is difficult to analytically calculate the numerator of 
(|24[1 . we resort to numerical integration. We can calcu- 
late also C{k) by the numerical integration of the prod- 
uct of C(o) and P(a\k) which is given by Bayes's law 
P(a\k) = P(k\a)p(a)/P{k) (see Fig. 2). For 7 > 2 the 
clustering coefficient C(k) decreases with the dimension 
d. On the other hand, for 7 < 2, the clustering coefficient 
C(k) seems to be independent of the dimension d. This 
suggests that when 7 < 2, the spatial structure is irrele- 
vant to the network structure. This fact is confirmed by 
the behavior of the average cluster coefficient defined by 
Watts and Strogatz 0, as is shown in Fig. 2 (d). 

Finally we study the average shortest path length L 
(Fig. 3). For 7 > 3, the average shortest path length 
seems to follow a power law L oc N^, where fi is some- 
what smaller than l/d. On the other hand, for 7 < 3, 
the average shortest path length grows more slowly than 
In TV. In this case, a pair of vertices with sufficiently 
large fitness are always linked, because max((9aiOj) — 
ON -1 ^> 1 for large N. Consequently some vertices, 
which connect each other regardless of their distance, 
compose a shortcut network. As a result, the spatial 
structure is irrelevant to the average shortest path length, 
and thus the network is ultrasmall [jj. 

In summary, we have studied a scale free network em- 
bedded on geographical space. While the scaling expo- 
nent of the degree distribution equals 7 for the scaling ex- 
ponent 7 of the fitness distribution for 7 > 2, it is always 
2 for 7 < 2. For 7 < 2, the spatial effect is irrelevant to 
some topological properties (ANND or clustering coeffi- 
cient) of the network. While the network is disassortative 
(negative degree correlation) for 7 < 3, it is weakly assor- 
tative (positive degree correlation) for 7 > 3. Moreover, 
the network is not small for 7 > 3, whereas the spatial 
effect is irrelevant to the average shortest path length for 
7 < 3. Thus, there are two crossovers at 7 = 2 and 7 = 3. 
The reason why the crossovers are observed clearly is 
that the fitness is assigned deterministically. If we use 
random fitness following the power law, the crossovers 
became somewhat blurred, but do not change qualita- 



tively. Furthermore, the preliminary numerical research 
suggests that these results hold even if the distribution of 
the location of the vertices is not uniform. Therefore, we 
expect that the results presented in this paper are robust 
in the condition that the nearer pairs tend to be linked. 

This research was carried out under the ISM Cooper- 
ative Research Program 2006-ISM CRP-1008. 
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